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ABSTRACT 

The rcdshifted spectral line radiation emitted from both atomic fine-structure and 
molecular rotational transitions in the interstellar medium (ISM) of high-redshift 
galaxies can be detected in the centimetre, millimetre and submillimetre wavebands. 
Here we predict the counts of galaxies detectable in an array of molecular and atomic 
lines. This calculation requires a reasonable knowledge of both the surface density of 
these galaxies on the sky, and the physical conditions in their ISM. The surface density 
is constrained using the results of submillimetre- wave continuum surveys. Follow-up 
OVRO Millimeter Array observations of two of the galaxies detected in the dust con- 
tinuum have provided direct measurements of CO rotational line emission at redshifts 
of 2.56 and 2.81. Based on these direct high-redshift observations and on models of 
the ISM that are constrained by observations of low-redshift ultraluminous infrared 
galaxies, we predict the surface density of line-emitting galaxies as a function of line 
flux density and observing frequency. We incorporate the sensitivities and mapping 
speeds of existing and future millimetre/submillimetre-wave telescopes and spectro- 
graphs, and so assess the prospects for blank-field surveys to detect this line emission 
from gas-rich high-redshift galaxies. 

Key words: ISM: molecules - galaxies: evolution - galaxies: formation - cosmology: 
observations ~ infrared: galaxies - radio lines: galaxies 



1 INTRODUCTION 

The redshifted far-infrared/submillimetre-wave line emis- 
sion from the interstellar medium (ISM) in galaxies could be 
exploited to detect new samples of distant gas-rich galaxies 
and active galactic nuclei (AGN) (Loeb 1993; Blain 1996; 
van der Werf & Israel 1996; Silk & Spaans 1997; Stark 
1997; Combes, Maoli & Omont 1999; van der Werf 1999). 
This emission is attributable both to molecular rotational 
transitions, in particular from carbon monoxide (CO), and 
to atomic fine-structure transitions, in particular from the 
singly ionized 158-pim carbon [Cll] line. 

Redshifted CO emission has been detected from a range 
of known high-redshift galaxies and quasars in the millime- 
tre/submillimetre waveband, as summarized by Frayer et al. 
(1998) and Combes et al. (1999). Many of these galaxies are 
known to be gravitationally lensed by a foreground galaxy, 
an effect which potentially complicates the interpretation of 
the results by altering the ratios of the inferred luminosities 



in the continuum and the detected lines (Eisenhardt et al. 
1996). In only a small subsample of the detected galaxies 
(Solomon, Downes & Radford 1992; Barvainis et al. 1997; 
Downes et al. 1999) have multiple lines been detected, pro- 
viding an opportunity to investigate the astrophysics of the 
ISM. 

So far there have been very few detections of redshifted 
fine-structure lines, despite careful searches, for example, 
for both [Cll] (Isaak et al. 1994; Ivison, Harrison & Coul- 
son 1998a; van der Werf 1999) and singly ionized 205-/im 
[Nil] emission (Ivison & Harrison 1996) . Neutral carbon [Cl] 
emission, which is considerably less intense than [Cll] and 
[Nil] emission in the Milky Way (Wright et al. 1991) and 
nearby galaxies (Stacey et al. 1991), has been detected from 
the gravitationally lensed Cloverleaf quasar (Barvainis et al. 
1994). The most luminous high-redshift galaxies and quasars 
have necessarily been targetted in these searches. 

[Cll] fine-structure emission is powerful in both the 
Milky Way and in sub-L* galaxies, in which it accounts for 
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about 0.5 per cent of the bolometric far-infrared luminosity 
(Nikola ct al. 1998). However, based on observations of a 
limited number of low-redshift galaxies using the Infrared 
Space Observatory (ISO) (Malhotra et al. 1997; Luhman et 
al. 1998; Piorini ct al. 1999), it appears that a systematically 
lesser fraction of the bolometric luminosity of more luminous 
galaxies appears as [Cii] emission, about 0.1 per cent (Luh- 
man et al. 1998). As noted by Luhman et al. (1998) and 
van der Werf (1999), the results of these ISO observations 
arc fully consistent with the non-dctcction of rcdshiftcd fine- 
structure emission from high-redshift gala^xies using ground- 
based submillimetre-wave telescopes. The results of Kuiper 
Airborne Observatory (KAO) observations of the Galactic 
centre (Erickson et al. 1991) indicate that 63- and 146-/im 
neutral oxygen [Ol] fine-structure emission becomes steadily 
more luminous as compared with that from [Cll] as the far- 
infrared luminosity of gas clouds increases. However, cur- 
rently there is insufficient published data available to ad- 
dress this issue in external galaxies. 

We attempt here to predict the counts of distant gas- 
rich line-emitting galaxies that could be detected in the mil- 
limetre/submillimetre waveband. There are two challenges 
to making reliable predictions. First, there are limited data 
available from which to construct a clear understanding of 
the astrophysics of the ISM in high-redshift galaxies. There 
are only a few tens of detections of line emission from these 
objects, the majority of which have been made in galaxies 
that arc gravitationally Icnscd by foreground galaxies. Be- 
cause of the potential for differential magnification across 
and within the lensed galaxy, neither the ratios of the line 
and continuum luminosities nor the excitation conditions in 
the ISM are known accurately in these cases. As shown by 
ISO [Cii] observations, extrapolation of the observed prop- 
erties of low-redshift galaxies with relatively low luminosi- 
ties to greater luminosities in high-redshift galaxies is not 
necessarily reliable. Secondly, the space density and form of 
evolution of gas-rich galaxies at high redshifts has not been 
well determined. Thus the existing predictions of the ob- 
servability of high-redshift submillimetre-wave line emission 
have concentrated on either discussing the potential observ- 
ability of individual high-redshift galaxies (van der Werf & 
Israel 1996; Silk & Spaans 1997; Combes ct al. 1999; van 
der Werf 1999), or have relied on extensive extrapolations, 
from the populations of low-redshift ultraluminous infrared 
galaxies (ULIRGs) (Blain 1996) and from low-redshift L* 
bulges to the properties of proto-quasars at ~ 10 (Loeb 
1993). 

Both of these difficulties can be addressed by exploiting 
the results of deep 450- and 850-/im dust continuum radi- 
ation surveys made using the Submillimetre Common-User 
Bolometer Array (SCUBA) camera (Holland ct al. 1999) at 
the James Clerk Maxwell Telescope (JCMT; Small, Ivison & 
Blain 1997; Barger et al. 1998; Hughes et al. 1998; Barger, 
Cowie & Sanders 1999; Blain et al. 1999b, 2000; Eales et 
al. 1999). These surveys are sensitive to galaxies at very 
high redshifts (Blain & Longair 1993), and have detected a 
considerable population of very luminous dust-enshrouded 
galaxies. The 15-arcsec angular resolution of the JCMT is 
rather coarse, but reliable identifications can be made by 
combining the SCUBA images with multi- waveband follow- 
up images and spectra (Ivison et al. 1998b, 2000; Smail et al. 
1998, 2000; Barger et al. 1999b; Lilly et al. 1999), and cru- 



cially with observations of rcdshiftcd CO emission, which are 
currently available for two submillimetre-selcctcd galaxies: 
SMMJ02399-0136 at z = 2.81 and SMM J14011-F0252 at 
z = 2.56 (Prayer et al. 1998, 1999; Ivison et al. 1998b, 2000). 
The bolometric and CO-line luminosities of these galaxies 
are rcEisonably well known, and because they are lensed by 
clusters rather than individual foreground galaxies, their in- 
ferred line ratios are not subject to modification by lens- 
ing. These observations thus provide a useful template with 
which to describe the properties of the ISM and line emission 
in high-redshift, dust-enshrouded, gas-rich galaxies. 

In Section 2 we discuss the existing line observations, 
and summarize our current state of knowledge about the 
evolution and rcdshift distribution of galaxies that have 
been discovered in submillimetre-wave dust continuum sur- 
veys. In Section 3 we describe our model of line emission 
from these galaxies, and present the results, as based on our 
understanding of high-redshift continuum sources. In Sec- 
tion 4 we discuss the observability of this hypothetical pop- 
ulation using existing and future millimctrc/submillimetre- 
wave spectrographs. Unless otherwise stated, we assume 
that Ho = 50kms"^ Mpc~^, fio = 1 and CIa = 0. 



2 BACKGROUND INFORMATION 

2.1 Line observations 

Ground-based telescopes have detected molecular rotation 
lines and atomic fine-structure linos from low-rcdshift galax- 
ies (e.g., Sanders et al. 1986; Wild et al. 1992; Devereux et 
al. 1994; Gerin & Phillips 1999; Mauersberger et al. 1999). 
Atomic fine-structure lines have also been observed from 
bright galactic star-forming regions and nearby galaxies us- 
ing the KAO (Stacey ct al. 1991; Nikola et al. 1998), COBE 
(Wright et al. 1991) and ISO (Malhotra et al. 1997; Luhman 
et al. 1998; Pierini et al. 1999). 

CO rotational line emission has been detected suc- 
cessfully from various high-redshift galaxies and quasars, 
including the first identified high-redshift ULIRG IRAS 
F10214-I-4724 (Solomon et al. 1992), the gravitationally 
lensed Cloverleaf quasar H 1413-1-117 (Barvainis et al. 1994; 
Kneib ct al. 1998), various quasars at 2: ~ 4, including 
BR 1202-0725 (Ohta et al. 1996, 1998; Omont et al. 1996) 
and the extremely luminous APM 08279-1-5255 (Lewis et al. 
1998; Downes et al. 1999), and the submillimetre-selected 
galaxies SMM J02399-0136 and SMM J14011-I-0252 (Prayer 
et al. 1998, 1999). A significant fraction of the dynamical 
mass in many of these systems is inferred to be in the form 
of molecular gas, and it is plausible that they are observed in 
the process of forming the bulk of their stellar populations. 

2.2 Atomic fine-structure lines 

Atomic fine-structure lines omitted at wavelengths longer 
than about 100 ^m - [Cll] at 1900 GHz, [Nil] at 1460 and 
2460 GHz, [Oi] at 2060 GHz, and [Ci] at 492 and 809 GHz 
- are redshifted into atmospheric windows for galaxies at 
redshifts z ^ 5, the rcdshift range within which at least 
80 per cent of dust-enshrouded galaxies detected by SCUBA 
appear to lie (Smail et al. 1998; Barger et al. 1999b; Lilly 
et al. 1999). There are many other mid- infrared lines with 
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shorter restframe emission wavelengths (see, e.g., Lutz et 
al. 1998); however, unless massive galaxies exist at z ~ 10, 
these lines will not be redshifted into atmospheric windows 
accessible to ground-based telescopes. 

Here we assume that the line-to-bolometric luminos- 
ity ratio /line = 10~* for the [Cll] line in all high-rcdshift 
gas-rich galaxies, corresponding to the value observed in 
low-redshift ULIRGs. The equivalent value for [Ci]492GHz, 
/lino = 2.9 X 10^® is chosen to match the value observed by 
Gerin & Phillips (1999) in Arp220. The values of /ino for 
other fine-structure lines listed in Table 1 are chosen by scal- 
ing the results of observations of the Milky Way and low- 
rcdshift galaxies (Gcnzcl ct al. 1990; Erickson ct al. 1991; 
Stacey et al. 1991; Wright et al. 1991). This approach should 
lead to a reliable estimate of the counts of [Cl] and [Cll] lines, 
but greater uncertainty in the N and O line predictions. In 
the absence of more observational data, which would ideally 
allow luminosity functions to be derived for each line, we 
stress that the predictions of the observability of redshifted 
fine-structure lines made in Sections must be regarded as 
tentative and preliminary. 



2.3 CO rotational transitions 

2.3.1 CO line excitation 

Much more observational data is available about the prop- 
erties of the ladder of CO rotational transitions in the ISM. 
The energy of the Jth level in the CO molecular rotation 
ladder Ej = k^T, where T = J( J + 1)[2.77K], and so the 
energy of a photon produced in the J + 1 ^ J transition 
is hv = fcB(J + 1)[5.54K]. The population of the J states 
can be calculated by assuming a temperature and density 
for the emitting gas. The primary source of excitation is ex- 
pected to be collisions with molecular hydrogen (H2), which 
dominates the mass of the ISM, with a role for radiative exci- 
tation, including that attributable to the cosmic microwave 
background radiation (CMBR). By taking into account the 
spontaneous emission rate, j4j+i,j oc {J / {2J +2i) with 
^1,0 = 6 X 10~^ s^^, and details of the optical depth and ge- 
ometry of gas and dust in the emitting region, the luminosi- 
ties of the various J -|- 1 — > J rotational transitions can be 
calculated. If the J state is to be thermally populated, then 
the rate of CO-H2 collisions in the ISM gas must be greater 
than about ^7+1 /• This condition will not generally be met 
for a temperature of 50 K in the CO (5^4) transition unless 
the density of H2 molecules exceeds about 2 x lO^cm"^, 
which is many times denser than the 10* cm~^ that appears 
to be typical of low-redshift ULIRGs (Downes & Solomon 
1998). Radiative excitation and optical depth effects, per- 
haps in very non-isotropic geometries, with very pronounced 
substructure, will complicate the situation greatly in real 
galaxies. In general, calculations of level populations are very 
complex, and at high rcdshifts there are very few data with 
which to constrain models. 

Probably the best way to investigate the conditions in 
very luminous distant galaxies is to study their low-redshift 
ULIRG counterparts, and the rare examples of high-redshift 
galaxies and quasars for which more than one CO transition 
has been detected (e.g. Downes et al. 1999), bearing in mind 
the potential eflfects of gravitational lensing. 

In order to try and make reasonable predictions for the 



line ratios in high-redshift galaxies, we employed a stan- 
dard large velocity gradient (LVG) analysis (e.g. de Jong, 
Dalgarno k, Chu 1975) to estimate how the CO line rar 
tios are affected by the temperature, density and finite 
spatial extent of the ISM, and by the radiative excita- 
tion caused by the CMBR. In the third column of Ta- 
ble 1 we show the results from this model, assuming a den- 
sity of lO^'cm"^, which is typical of the central regions of 
ULIRGs (Downes & Solomon 1998), and a standard value 
of X(CO)/(d'(j/dr) = 3 x lO"'' (kms-ipc-i)-^ We as- 
sume a kinetic temperature of 53 K, which is the temper- 
ature of the dominant cool dust component in the SCUBA 
galaxy SMM J02399-0136 (Ivison ct al. 1998b). Higher dust 
temperatures of about 80 and 110 K are inferred for other 
well-studied high-redshift galaxies IRAS F10214-I-4724 and 
APM 08279-1-5255 respectively, but these are very exotic 
galaxies, and the results are potentially modified by the ef- 
fects of differential gravitational lensing. We assume a back- 
ground temperature of 10 K, the temperature of the CMBR 
at z — 2.7, the mean redshift of the two SCUBA galaxies 
with CO detections. The line-to-continuum bolometric lu- 
minosity ratio /line = Lj+i^j/Lpin can be calculated if the 
bolometric continuum luminosity Lfir is known. There is a 
clear trend of a reduction in the CO line-to-bolometric lumi- 
nosity ratio of luminous infrared galaxies as the bolometric 
luminosity increases, with a large scatter, which is consis- 
tent with /i 

inc OC ^pjpj^ (Sanders et al. 1986). AVe normalize 
the results to the observed CO (3^2) line luminosity in the 
LpiR ~ 10" Lq scuba galaxies SMM J02399-0136 and 
SMM J14011-I-0252, in which the ratios of the luminosity in 
the CO(3^2) line to Lfir is about 2.1 x 10"® and 5.3 x 10~® 
respectively (Fraycr et al. 1998, 1999), with errors of order 
50 per cent. The dominant source of error is the uncertainty 
in the bolometric luminosity. 

We compare the results obtained in a simple equilib- 
rium case, in which the density of CO molecules, the spon- 
taneous emission rates and the transition energies in differ- 
ent J states are multiplied to give the luminosity in each 
transition, 

Lj+i^j oc i/^(J + if exp{-[2.77K](J + 1)(J + 2)/T}. (1) 

The results are shown in the fourth and fifth columns of 
Table 1 for J < 9, assuming kinetic temperatures of 38 and 
53 K respectively, the temperature generated by simple fits 
to the observed counts of IRAS and ISO galaxies (Blain et 
al. 1999c), and the spectral energy distribution (SED) of 
SMM J02399-0136. 

2.3.2 The effects of different excitation conditions 

The values of func listed in Table 1 for CO transitions in the 
LVG and the 38- and 53-K thermal equilibrium models dif- 
fer. However, in lines with J g 7 the differences are less than 
a factor of a few. Given the current level of uncertainty in the 
data that support these calculations, this is an acceptable 
level. The differences between the results are more marked 
at large values of J. The consequences of these differences 
for the key predictions of the source counts of line-emitting 
galaxies arc discussed in Section 3.2. 

Throughout the paper we use the LVG model to 
describe the CO line emission of dusty galaxies. Ob- 
servations of the CO line ratios in low-redshift dusty 
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Table 1. The fraction of the bolometric luminosity of distant dusty galaxies that is assumed to be emitted in a variety of submillimetre- 
wave lines /unei their restframe emission frequencies i^rest, and the redshifts at which the lines would be detected in the important 
observing bands at 90, 230, 345 and 650 GHz; zgo, Z230, ^345 and zeso respectively. The CO ratios /une listed in the third, fourth and 
fifth columns are calculated assuming the LVG model described in Section 2, and in two local thermal equilibrium models with kinetic 

temperatures of 38 and 53 K respectively. The fraction of the bolometric luminosity of the galaxies in all CO transitions, obtained by 
adding all values of fn^c ^1^° shown. The values of fu^e listed for fine-structure transitions in the final six rows are derived from 
observations of low-rcdsliift galaxies and the Milky Way: see Section 2.2. The line styles and thicknesses used to represent the various 
transitions are listed in the final column. 
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galaxies indicate a wide range of excitation conditions, 
Tb[CO(3^2)/CO(1^0)] ~ 0.2-1 (Mauersberger et al. 
1999). In the central regions of starburst nuclei this tem- 
perature ratio tends to be systematically higher, with 
Tb[CO(3-»2)/CO(l-»0)] ~ 0.5-1 (Devereux et al. 1994). In 
our chosen LVG model, listed in Table 1, this ratio is ~ 0.9, 
which is consistent with the observations of the central re- 
gions of M82 (Wild ot al. 1992) and Arp 220 (Mauersberger 
et al. 1999). The value of the X parameter in the model has 
little effect on these ratios; however, reducing the density 
from 10'' to 3300 and 1000 cm reduces the predicted ratio 
to 0.81 and 0.51 respectively. Our model looks reasonable in 
the light of these observations, as the liigh-redshift galaxies 
would typically be expected to be ULIRGs with high gas 
densities. 

There have been two recent discussions of the observ- 
ability of CO line emission from high-redshift galaxies. Silk 
& Spaans (1997) describe the effect of the increasing radia- 
tive excitation of high- J lines at very high redshifts because 
of the increasing temperature of the CMBR. The median 
redshift of the SCUBA galaxies is likely to be about 2-3 
(Bargcr ct al. 1999b; Small et al. 2000; Lilly et al. 1999), 
with perhaps 10-20 per cent at z ^ 5, and because there is 
currently no strong evidence for the existence of a large pop- 
ulation of metal-rich galaxies at 2; ^ 10, this effect is unlikely 
to be very important. Combes et al. (1999) include a hot 
dense 90-K/lO'' cm-'^ component in the ISM of their model 
high-redshift galaxies, in addition to the cooler less dense 
component included in our models, and conduct LVG calcu- 
lations to determine CO emission-line luminosities. Under- 



standably, the luminosity of high- J CO lines is predicted to 
be greater in their models as compared with the values listed 
in Table 1. The continuum SED of the best studied SCUBA 
galaxy SMM J02399— 0136 certainly includes a contribution 
from dust at temperatures greater than 53 K, but here we 
avoid including additional hot dense phases of the ISM in 
our models in order to avoid complicating the models and to 
try and make conservative predictions for the observability 
of high- J CO lines at high redshift. 

Only additional observations of CO in high-redshift 
galaxies will allow us to improve the accuracy of the con- 
ditions in the ISM that are assumed in these models. The 
detection of the relative intensities of the C0(9— >8) and 
CO(5^4) emission from APM 08279+5255 (Downes et al. 
1999), and the ratio of the intensities of the multiple CO 
lines detected in BR 1202-0725 at z = 4.7 (Ohta ct al. 1996, 
1998; Omont et al. 1996) are broadly consistent with the val- 
ues of /line Usted in the 53-K thermal equilibrium model: see 
Table 1. 

The evolution of the abundance of CO and dust 
throughout an episode of star formation activity have been 
investigated by Fraycr & Brown (1997), and the detailed 
appearance of the submillimetre-wave emission-line spec- 
trum of galaxies requires a careful treatment of the radiative 
transfer between stars, AGN, gas and dust in an appropriate 
geometry. However, given that the amount of information 
on the spectra of high-redshift galaxies is currently not very 
great, it seems sensible to base estimates of the properties of 
line emission on the template of the submillimetre-selected 
galaxies studied by Prayer et al. (1998, 1999). 
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2.3.3 Other molecular emission lines 

There could also bo a contribution from rotational lines 
emitted by other species, such as NH3, CS, HON, HCO'*' 
and H2O; however, it seems unlikely that these emission 
lines would dominate the energy emitted in CO unless the 
densities and excitation temperatures are very high. 



2.4 High-redshift dusty galaxies 

The surface density of 850-/im SCUBA galaxies is now 
known reasonably accurately between flux densities of 1 and 
lOmJy (Barger et al. 1999a; Blain et al. 1999b, 2000). By 
combining knowledge of the properties of the SCUBA galax- 
ies, the low-redshift 60-/Lim IRAS galaxies (Saunders et al. 
1990; Soifcr & Ncugobaucr 1991), the 175-^im ISO galax- 
ies (Kawara et al. 1998; Puget et al. 1999) and the inten- 
sity of far-infrared background radiation (Fixsen et al. 1998; 
Hauser et al. 1998; Schlegel, Finkbeiner & Davis 1998), it is 
possible to construct a series of self-consistent models that 
can account for all these data (Guiderdoni et al. 1998; Blain 
et al. 1999b, c), under various Eissumptions about the forma- 
tion and evolution of galaxies. The 'Gaussian' model, which 
is described by Blain et al. (1999b) and based on pure lumi- 
nosity evolution of the low-redshift 60- /im luminosity func- 
tion (Saunders et al. 1990), was modified slightly to take 
account of the tentative redshift distribution derived from 
the SCUBA lens survey by Barger et al. (1999b). This 'mod- 
ified Gaussian' model is used as a base for the predictions 
of line observability presented here. The evolution function 
increases as (I-I-2;)''' with 7 ~ 4 at z ;S 1, has a 1.0-Gyr-long 
Gaussian burst of luminous activity centred aX z = 1.7, in 
which L* is 70 times greater than the value of L* at « = 
(Saunders et al. 1990), and then declines at 2; ^ 3. 

This model contains fewer parameters than there arc 
constraining pieces of information, and so should provide 
a reasonable description of the properties of high-redshift 
dust-enshrouded galaxies, which are expected to be the most 
easily detectable sources of line omission. Future observa- 
tions will inevitably provide more information and demand 
modifications to the model of galaxy evolution discussed 
above; however, it is likely to provide a sound basis for the 
predictions below. 



3 LINE PREDICTIONS 

In this paper we are concerned with the detcctability of red- 
shifted linos rather than with their resolution. As a result, we 
want to estimate the total luminosity of a line, and will not 
be concerned with details of its profile. The results are all 
presented as integrated flux densities, determined over the 
whole line profile. Where relevant, a line width of 300 km s~^ 
is assumed. 

The evolution of the 60- /im luminosity function of 
dusty galaxies (Saunders et al. 1990) is assumed to be de- 
fined by the modified Gaussian model. The bolometric far- 
infrared continuum luminosity function $boi(iFiR ) can then 
be calculated by integrating over a template dusty galgixy 
SED (Blain et al. 1999b). $boi can in turn be converted 
into a luminosity function for each line listed in Table 1, 



*l?iinG(iiino, z), by evaluating <l>boi at the bolometric lumi- 
nosity Lfir that corresponds to the line luminosity Lune. 

Based on observations by Sanders et al. (1986), the ra- 
tio of the luminosity in the CO (1^0) line, Leo, to Lfir is a 
function of Lfir, with Leo <x Lpf^. Hence, when making the 
transformation from Lune to Lfir for the CO lines listed in 
Table 1, we use the relationship Lfir = Lune/Zune, in which 
/line oc Lp°^ and is normalized to the value listed in Ta^ 
ble 1 at 10^^ Lq, the luminosity of the submillimetre-selected 
galaxies SMM ,102399-0136 and SMM J14011+0252. There 
is no evidence from ISO observations of any systematic lu- 
minosity dependence in the value of /une for [Cii] emission 
from ULIRGs, and so we assume that the values of func 
listed in Table 1 describe the linc-to-bolomctric luminosity 
relation in the fine-structure lines at all luminosities, that is 
Lfir = Liine//iine, where /une is constant. 

The surface density of line-emitting galaxies N{> S) 
that can be detected at integrated flux densities brighter 
than S, as measured in Jykms^^ or Wm^'^, in an observ- 
ing band spanning the frequency range between Vohs and 
Vobs + Aj/obs can be calculated by integrating the luminosity 
function of a line $iine over the redshifts for which the line 
is in the observing band, and luminosities Lunc greater than 
the detection limit Lmin(>S', «)• The count of galaxies is thus 

f.Z2 roc , 
N{>S)= / / $li„e(Lline,«)dLli„eD'(,?)^dz. (2) 

D is the comoving distance parameter to redshift z, and dr 
is the comoving distance element. For a line with a restframe 
emission frequency I'une, the limits in equation (2) are 

z\ = 1, (3) 

Inline 

and 

Z2 = 1. (4) 

t'linc 

If Zi > «o, where Zq is the maximum redshift assumed for 
the galaxy population, or if Z2 < 0, then the count N = 
0. Here zq — 10 is assumed. The minimum detectable line 
luminosity, 

Lrr^US, Z) = 477S'L»'(1 + zf . (5) 

If the integrated flux density of a line S observed at fre- 
quency V, is lWm~^, then the same quantity can be ex- 
pressed as 3(j//Hz)"^ X 10^^ Jykms"^ 

The counts calculated for the transitions listed in Ta- 
ble 1 are shown in Figs 1-3. Predictions are made at cen- 
tre frequencies (bandwidths) of 23 (8) - the radio K band, 
50 (8) - the radio Q band, 90 (1), 230 (8), 345 (8), 650 
(8), 200 (100) GHz, and 1 (l)THz, in Figs 2(a), 2(b), 1(a), 
1(b), 1(c), 1(d), 3(a) and 3(b) respectively. The 8-GHz radio 
bandwidth is matched to the performance specified for the 
upgraded VLA. The 1-GHz bandwidth at 90 GHz is approx- 
imately matched to the current performance of millimetre- 
wave interferometer arrays. The 8-GHz bandwidth at 230, 
345 and 650 GHz matches that for the SPIFI Fabry - Perot 
spectrograph (Stacey et al. 1996; Bradford et al. 1999, in 
preparation) and is a plausible value for the bandwidth of 
the future ground-based Atacama Large Millimetre Array 
(ALMA; Ishiguro et al. 1994; Brown et al. 1996; Downes et 
al. 1996), although the current goal is for a 16-GHz band- 
width (Wootten 2000). The count predictions in the very 
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Figure 1. Predicted counts of CO rotation and fine-structure lines detectable from galaxies in a 1-GHz band centred at 90 GHz 

(a) , and in 8-GHz bands centred on (b) 230 GHz, (c) 345 GHz and (d) 650 GHz. The band in (a) is accessible to the existing BIMA, 
IRAM, Nobeyama and OVRO interferometer arrays. These arrays also operate with a narrower 1-GHz bandwidth around 230-GHz 

(b) . ALMA will operate in all the bands shown in (b), (c) and (d), with a current goal bandwidth of 16 GHz (Wootten 2000). 
SPIFI and the HIFI and SPIRE-FTS instruments on the FIRST satellite will operate in the band shown in (d). The line styles and 
thicknesses correspond to each transition listed in Tablet. In (a) the CO rotational transitions from CO(l^O) to CO(8— >7) and 
the [Cl]492GHz line arc present, in (b) CO(2^1) to CO(9^8) and all the fine-structure lines listed in Tablet arc detectable, in (c) 
CO(3^2) to CO(9 ^8) and all the fine-structure lines are present, and in (d) the CO transitions from CO(6^5) to CO(9^8) and 
all but the [Ci]492Giiz fine-structure line are present. 10~^''Wm~^ is equivalent to 3.3, 1.3, 0.9 and 0.5 Jykms"^ at 90, 230, 345 
and 650 GHz respectively. 



wide atmospheric window between 150 to 250 GHz is shown 
in Fig. 3(a). This band cannot be observed simuhaneously 
using heterodyne instruments, but could plausibly be cov- 
ered using an advanced grating or Fabry-Perot spectrograph 
feeding a sensitive bolometer detector array. In Fig. 3(b) the 
far-infrared/submillimetre band between 460 and 1500 GHz 
is shown. This band is matched to the specified spectral 
range of the SPIRE Fourier Transform Spectrograph des- 



tined for the FIRST satellite (Griffin 1997; Griffin et al. 
1998). 



3.1 The shape of the counts and optimal surveys 

The integrated counts, N{> S) oc S~" presented in Figs 
1 to 3 all have a characteristic form, with a relatively flat 
slope, a ~ 0.3, at faint flux densities, which rolls over to 
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Figure 2. The line counts for galaxies at z < 10 predicted in a 8-GHz band centred at (a) 23 GHz and (b) 50 GHz, corresponding to 
the K and Q bands. In (a) only the CO{1^0) and CO(2^1) lines arc rcdshiftcd into the band. In (b) CO lines up to CO(5— >4) and 
the [Ci]492 GHz line are present. The line styles and thicknesses that correspond to the transitions are listed in Table 1. 10~^^ Wm~-^ 
is equivalent to 0.13 and 0.06 Jykms"^ at 23 and 50 GHz respectively. 

(a) (b) 




Integrated line flux density / ¥ m ^ Integrated line flux density / W m" 



Figure 3. The line counts predicted in very wide bands covering the spectral ranges between 150-250 GHz (a) and 460-1500 GHz 
(b). The band in (a) could be observed using a grating spectrograph feeding bolometer detectors. The band in (b) could be 
surveyed using the SPIRE-FTS instrument fitted to FIRST. The lowest CO transition shown in (a) is CO(2— *1). In (b) the lowest 
CO transition is CO(5— »4), and the counts are dominated by fine-structure lines. The line styles and thicknesses that correspond 
to each transition are listed in Table 1. 10"^" Wm~^ is equivalent to 1.5 and 0.4Jykms~^ at 200 and 800GHz respectively. 



a steep decline, a ~ 3, at brighter fluxes. Note that the 
detection rate of galaxies in a survey is maximized at the 
depth at which a = 2. Because the counts of lines presented 
here have a pronounced knee at a certain flux density, at 
which the value of a crosses 2, this flux density is the opti- 
mal depth for a blank-field line survey, and surveys that are 
either shallower or deeper should be much less efficient. 

In Fig. 4 the predicted counts of line-emitting galaxies 
are shown as a function of observing frequency between 50 



and 2000 GHz, at both faint and relatively bright values of 
the integrated line flux density S, 10"^^ and 10"^" Wm~^. 
The form of the curve for each emission line is determined 
by the form of evolution of fax-infrared luminous galaxies 
specifled in the modifled Gaussian model. In any particu- 
lar line, galaxies detected at higher frequencies are at lower 
redshifts. This effect accounts for the double-peaked distri- 
bution that is most apparent in the counts of bright galaxies 
shown in Fig. 4(b). The broad peak at lower frequencies cor- 
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(a) (b) 



T 




Frequency / GHz Frequency / GHz 

Figure 4. The surface density of line emitting galaxies brighter than (a) 10~^^ Wm~^ and (b) 10"^" Wm~^, as a function of 
observing frequency, assuming an observing bandwidth of 8 GHz. The line styles are listed in Table 1. The sharp lower cutoff to 
the counts of fine-structure lines at low frequencies is attributable to the assumed upper limit to the redshift distribution of dusty 
galaxies, zo = 10. If gas-rich galaxies exist at z > 10, then they would be detectable at frequencies below these cutoffs. 



responds to the very luminous Gaussian burst of activity at 
z ~ 1.7, while less luminous low-redsliift galaxies account 
for the second sharper peak at higher frequencies. The dou- 
ble peak is more noticeable for the brighter counts, to which 
the contribution of low- redshift galaxies is more significant. 

Based on the predicted counts of CO lines, which are 
probably quite reliable, the most promising frequency range 
in which to aim for CO line detections is 200-300 GHz, re- 
gardless of whether faint or bright CO lines are being sought. 
At frequencies greater than about 500 GHz the counts are 
expected to be dominated by fine-structure lines, and al- 
though these counts are still very uncertain, the dctcctabil- 
ity of fine-structure line emitting galaxies is expected to be 
maximized at a frequency between about 400 and 800 GHz. 
Future blank-field searches for fine-structure emission from 
distant galaxies should probably be concentrated in this 
range. 

3.2 The effects of CO excitation 

The three different sets of values of /une listed for CO tran- 
sitions in Table 1 were derived assuming different excitation 
conditions in the ISM. As the detailed astrophysics of gas- 
rich high-redshift galaxies is very uncertain and only weakly 
constrained by observations, there is an inevitable uncer- 
tainty in the counts predicted in Figs 1-4. 

The most significant differences between the values of 
/lino predicted in the LVG, 38-K and 53-K models occur in 
the higher J lines, for which the effect of the higher excita- 
tion temperature and the subthermal excitation of lines in 
the LVG model is greater. The increase is most significant at 
bright flux densities. If the alternative values of /une listed 
in the thermal equilibrium models arc used to predict counts 
similar to those in Figs 1-4, then for J < 7, at interesting line 
flux densities of about 10^^^ Wm~^, the counts differ by less 
than a factor of about 3. This uncertainty is greater than 



that in the surface density of high-redshift submillimetre- 
selected galaxies that were used to normalize the model of 
galaxy evolution, and so the counts predicted in Figs 1-4 
should be reliable to within a factor of about 3. Note that 
the counts in Fig. 4 indicate that the C0(5— >4) line is likely 
to be provide the most significant contribution to the counts, 
and so the degree of excitation of lines with J > 7 is not 
likely to affect the detection rate in future blank-field line 
surveys. 



3.3 The effects of cosmology 

In Fig. 5 the 230-GHz counts expected in three different 
world models are compared, based on the same underlying 
population of distant dusty galaxies. In world models with 
Qo = 0.3 and either Qa = 0.0 or 0.7 the faint counts are ex- 
pected to be greater, and the bright counts are expected to 
be less, as compared with the predictions of an Einstein-de 
Sitter model. The count changes as a function of cosmology 
by a factor of about 10 in the most abundant lines over the 
wide range of flux densities presented in Fig. 5, which cor- 
respond to a range of count values between about 10* and 
10"^ deg"^ At a flux density of 10"^° Wm'^ which is likely 
to be similar to the limiting depth of future submillimetre- 
wave line surveys, the counts are modified by only a factor of 
about 2 by changing the world model. Hence, at present the 
uncertainties in the excitation conditions in the ISM are ex- 
pected to be greater than those introduced by uncertainties 
in the cosmological parameters. 



4 LINE OBSERVATIONS 
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Figure 5. The effects of different world models on the predictions of CO and fine-structure line counts, for demonstration in a 

narrow band centred on 230 GHz (sec Fig. 2). The line styles and thicknesses that correspond to each transition are listed in Table 1. 
The results arc presented as ratios between the counts in two world models with f2o < 1, and the predictions of an Einstein— de 
Sitter model. In (a) the cosmological parameters arc SIq = 0.?> and Ha = 0.0; in (b) Qq = 0.3 and Q/^ = 0.7. The differences 
due to the world model are by factors of a few, at worst comparable with the level of uncertainty in the count predictions that is 
attributable to the uncertain excitation conditions in the ISM. 



4.1 Millimetre-wave interferometer arrays 

Four millimetre- wave interferometer arrays are currently op- 
erating: the BIMA array with ten 6-m antennas; the IRAM 
array with five 15-m antennas; the Nobeyama Millimeter 
Array with six 10-m antennas; and the OVRO Millimeter 
Array with six 10.4-m antennas. These instruments operate 
at wavelengths longer than 1 mm, and their fields of view and 
sensitivities translate into similar mapping speeds to equiva- 
lent flux density levels. For example, in a 20- hr integration at 
a frequency of 90 GHz, the OVRO Millimeter Array is able 
to reach a 5ct sensitivity limit of about 5 x 10~^^Wm~^ 
across a 1-GHz band in a 1-arcmin^ primary beam. 

The counts of lines expected in a 90-GHz observation 
with a 1-GHz bandwidth are shown in Fig. 1(a) as a function 
of integrated line flux density. At this limiting flux density, a 
count of about 2 deg~^ is expected, and so the serendipitous 
detection of an emission line in a blank-field survey would 
be expected to occur after approximately 4 x 10* hr. In 230- 
GHz observations with the same 1-GHz bandwidth, which 
are possible in good weather, the sensitivity required for a 5a 
detection in a 20-h integration is about 2 x 10~^''Wm~^, 
in a 0.25-arcmin'^ primary beam. The 230-GHz count ex- 
pected at this depth, after correcting the results shown in 
Fig. 1(b) for the narrower bandwidth, is about 50deg~^. A 
serendipitous line detection would thus be expected about 
every 6000 hr. Thus, while known high-redshift galaxies can 
certainly be reliably detected in reasonable integration times 
using these arrays, blank-field surveys for emission lines are 
not currently practical. 

The detection of rodshiftod CO lines at a 90-GHz 
fiux density of about 10~^^ Wm~^ from two submillimetre- 
wave continuum sources, with a surface density of several 
100 deg"^ (Blain et al. 1999b, 2000), using the OVRO ar- 



ray (Frayer et al. 1998, 1999) is consistent with the 90-GHz 
line count of about 2 deg~^ predicted in Fig. 1(a) at this fiux 
density. This is because observations in many tens of 1-GHz 
bands would have been required to search for a CO emission 
lino from those galaxies without an optical rodshift, and so 
the source count of lines at the detected fiux density is ex- 
pected to be about two orders of magnitude less than that 
of the dust continuum-selected galaxies. 

The large ground-based millimetre/submillimetre-wave 
interferometer array ALMA will provide excellent subarcsec- 
ond angular resolution and a large collecting area for obser- 
vations of submillimetre-wave continuum and line radiation. 
Based on the performance described for ALMA at 230 GHz 
by Wootten (2000), a 300-km-s~^ line with an integrated 
fiux density of 5 X 10^^^ Wm^'^ could be detected at 5a 
significance, but not resolved, anywhere within a 16-GHz 
band in about Ihr in the 0.15-arcmin^ primary beam. The 
surface density of lines brighter than this flux density is ex- 
pected to be about 2.5 x 10* deg^'^ (Fig. lb; see also Blain 
2000), and so a detection rate of about 1.7 hr~^ would be 
expected. The knee in the counts in Fig. 1(b), indicating the 
most efficient survey depth, is at a flux density greater than 
about 10-2" W m ^, a depth reached in an integration time 
of about 10 s at 5a significance. Hence, making a largo mo- 
saic map at this depth, covering an area of 0.017 dog^ hr"^ 
should maximize the detection rate, which would then be 
about 15 galaxies per hour, neglecting scanning overheads 
incurred in the mosaiking process. This detection rate would 
allow large samples of line-emitting high-redshift galaxies 
to be compiled rapidly using ALMA. The performance of 
ALMA in line searches is discussed in more detail by Blain 
(2000). 
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4.2 Ground-based single-antenna telescopes 

At present the heterodyne spectrographs fitted to the 
JCMT, the Caltech Submillimeter Observatory (CSO), and 
the IRAM 30-m antenna are not sufficiently sensitive to al- 
low a blank-field survey to search for distant line emitting 
galaxies. For example, the 230-GHz receiver at the JCMT, 
which is the least susceptible to atmospheric noise, can reach 
a 5(7 sensitivity of about 2 x 10~^^ Wm^'^ in a 1-hr obser- 
vation in a 1.8-GHz band centred on 230 GHz. At this flux 
density the surface density of CO line emitting galaxies is 
expected to be ^ lOdeg"'^ (see Fig. lb), and so because the 
beam area is 2.4 x 10~^ deg^, many tens of thousands of 
hours of observation would be required to detect a source 
serendipitously. The development of wide-band correlators, 
such as the 3.25-GHz WASP (Isaak, Harris & Zmuidzinas 
1999), will improve the performance of single-antenna tele- 
scopes significantly for the detection of faint lines in gala:xies 
with a known redshift, but unless bandwidths are increased 
by at least an order of magnitude, blank-field line searches 
from the ground will not be practical. 

4-2.1 Future instrumentation 

At present, submillimetre-wave emission lines from a partic- 
ular distant galaxy can only be detected if an accurate red- 
shift has been determined, because the instantaneous band- 
width of receivers is narrow, and accurate tuning is required 
for a target line to be observed within the available band. For 
galaxies and AGN detected in blank-field submillimetre con- 
tinuum surveys, obtaining a spectroscopic redshift requires 
a very considerable investment of observing time (see, e.g.. 
Small et al. 1999a). However, if a rrmch wider bandwidth 
of order 115 GHz/ (1 -|- z) could be observed simultaneously, 
then a CO line would always lie within the observing band 
from a galaxy at redshift z. 

There are currently two ways to increase the instan- 
taneous bandwidth of a ground-based millimetre-wave tele- 
scope, both based on bolometer detectors rather than het- 
erodyne mixers. Either a resonant cavity, such as a Fabry- 
Perot, or a diffraction grating could be used to feed an array 
of bolometer detectors. The potential of such instruments 
are discussed in the subsections below. 

4.2.2 A Fabry-Perot device: SPIFI 

SPIFI is a 5 X 5 element Fabry-Perot interferometer, for 
use at frequencies from 460 to 1500 GHz on both the 1.7-m 
AST-RO telescope at the South Pole and the 15-m JCMT 
(Stacey et al. 1996; Bradford ct al. 1999, in preparation). On 
AST-RO, in the 350-/im (850-GHz) atmospheric window, 
the field of view of the instrument is about 25 arcmin^ . At a 
coarse 300-km-s~^ resolution, a bandwidth of about 8 GHz 
can be observed to a 5a detection threshold of about 3.2 x 
10"^^ Wm"^ in a 1-hr integration. On the 15-m JCMT, at 
the same frequency and for the same bandwidth, the field 
of view and the 5a detection threshold are both less, about 
0.3arcmin^ and 8 x 10~^® Wm~^ respectively. 

From Fig. 1(d), the count at flux densities brighter than 
the 5a detection threshold for the JCMT are expected to 
be about 30deg~^, and to be dominated by [Cii]-emitting 
galaxies, for which the counts are very uncertain. Hence 



a serendipitous detection would be expected about every 
400 hr, and so blank-fleld line surveys using SPIFI are on 
the threshold of being practical. A second-generation instru- 
ment with a wider field of view should be capable of making 
many detections in a reasonable integration time. 

4.2.3 A millimetre-wave grating spectrograph 

An alternative approach to obtaining a wide simultane- 
ous bandwidth would be to build a grating spectrograph 
to disperse the signal from a 10-m ground-based tele- 
scope onto linear arrays of bolometers. Such an instrument 
would be able to observe a reasonably large fraction of the 
clear 150-250 GHz atmospheric window simultaneously, at 
a background-limited la sensitivity of about 10^^° Wm~^ 
in a 1-hr observation, which is uniform to within a factor of 
about two across this spectral range. The predicted counts 
of lines in this spectral range are of order 100 deg~^ at this 
depth (Fig. 3a), and are expected to be dominated by CO 
lines, for which the predicted counts should be reasonably 
accurate. At 200 GHz the flcld of view of a 10-m telescope 
is about 0.2arcmin^, and so a 5a detection rate of about 
0.05 hr~^ would be expected in a blank-field line survey. A 
grating instrument with these specifications is thus on the 
threshold of being useful for conducting such a survey. 

It would also be very valuable for detecting CO lines 
emitted by the high-redshift galaxies detected in dust con- 
tinuum surveys whose positions are known to within about 
10 arcsec. An instantaneous bandwidth of order 100 GHz 
would accommodate either 2 or 3 CO lines for a, z > 2 galaxy, 
and so a spectroscopic redshift could be determined with- 
out recourse to radio, optical or near-infrared telescopes. For 
example, the 2 = 2.8 SCUBA galaxy SMM J02399-0136 
has an integrated flux density of 1.5 x 10~^^Wm~^ in 
the CO (3^2) transition (Frayer et al. 1998), and all the 
CO( J-l-1— > J) transitions with J=5, 6, 7 and 8 fall within the 
150-250 GHz spectral range. Based on the LVG model (Tar 
ble 1), these lines are expected to have integrated flux densi- 
ties of 3.0, 3.6, 3.2 and 1.0 x 10"^" W m"^ respectively. Thus, 
in a 5-hr integration in the fleld of SMM J02399-0136, aU 
but the C0(8— >7) line could be detected using a background- 
limited grating spectrograph, allowing its redshift to be de- 
termined unequivocally in a comparable time to that re- 
quired to obtain an optical redshift using a 4-m telescope 
(Ivison et al. 1998b). There axe two further advantages. 
First, the CO-line redshift would be the redshift of cool gas 
in the ISM of the galaxy, and not that of optical emission 
lines, which are typically blueshifted by several lOOkms"^ 
with respect to the ISM. Secondly, the ratios of the luminosi- 
ties of the detected CO lines would reveal information about 
the physical conditions in the cool ISM where star-formation 
takes place. 

4.3 Ground-based centimetre-wave telescopes 

Low- J CO lines are redshifted into the centimetre waveband 
for redshifts z < 10. Instruments operating in this waveband 
include the VLA and the Green Bank Telescope (GBT). 

Currently, the VLA can observe spectral lines in a very 
narrow 87.5-MHz band, centred in the K band (22-24 GHz) 
and Q band (40-50 GHz). In the K band, Ivison et al. (1996) 
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attempted to detect CO(l^O) emission from the environ- 
ment of a dusty radio galaxy at z = 3.8, and in 12 hours 
they obtained a 5a upper limit of about 4 x 10~^^ Wm~^ in 
a 2-arcmin beam. Because of the very narrow bandwidth, the 
detection rate of CO line-emitting galaxies with the K-band 
VLA in a blank- field survey is expected to be about lyr~^. 
However, by 2002 fibre-optic Unks and a new correlator will 
be installed, increasing the simultaneous bandwidth greatly 
to 8 GHz, and improving the 5(7 sensitivities to 1.3 x 10~^^ 
and 8.0 x 10^'^'^ Wm^'^ in 12-hr K- and Q-band integrations 
respectively. The counts of line-emitting galaxies at « < 10 
in these bands are shown in Fig. 2. At these sensitivities, 
K- and Q-band counts of about 500 and 2000 dcg"'^ arc ex- 
pected respectively, each corresponding to a detection rate 
of about 0.03 hr~^. The knee in the counts at which the de- 
tection rate is maximized is expected at rather similar flux 
densities of about 3 x 10~^^ and 2 x 10~^^ Wm~^ in the 
K and Q bands respectively, depths which can be reached 
in about 140- and 190-rriin integrations and at which the 
counts are expected to be about 100 and 400dcg^~ respec- 
tively. The most efficient detection rates in the K and Q 
bands are thus expected to be about 0.04 and 0.03 hr~^. 
Hence, low- J blank- field CO-line surveys will be practical 
using the upgraded VLA. High-redshift galaxies with lu- 
minosities comparable to the SCUBA galaxies detected by 
Prayer et al. (1998, 1999) should be readily detectable, as 
their integrated flux densities in the CO(l^O) line would 
be of order 7 x 10"^'^ and 4 x 10"^^ Wm"^ when redshifted 
into the K and Q bands, from z = 4 and 1.3 respectively. 

The 100-m clear-aperture GBT will operate with a 
3.2-GHz bandwidth in the K and Q bands, reaching 5cr 
sensitivity limits in 1-hr integrations of about 1.3 and 
5 X 10^^^ Wm^^ respectively. These sensitivities make the 
GBT even more suitable for detecting low-J CO lines from 
known high-redshift galaxies than the VLA, but the subar- 
cminutc field of view of the GBT is too small to allow a 
practical blank-field survey. At these sensitivity limits, only 
2 X 10^^ and 6 x 10~* lines per beam are expected in the K 
and Q bands respectively. 

4.4 Air- and space-borne instruments 

The 3.5-m space-borne telescope FIRST {Pilhratt 1997) will 
carry the HIFI far- infrared/submillimetre- wave spectrome- 
ter (Whyborn 1997) and the SPIRE bolometer array cam- 
era (Griffin et al. 1998). Both of these instruments operate 
at frequencies for which fine-structure lines are expected to 
dominate the counts of line-emitting galaxies, and so the 
number of detectable galaxies is necessarily uncertain. 

The spectral coverage of HIFI extends from 500 to 
1100 GHz. The instrument will have a 4-GHz bandwidth, 
a 0.4-arcmin^ field of view and a 5a sensitivity of about 
5 X 10"^^Wm~^ in a 1-hr integration at 650 GHz. The 
SPIRE Fourier Transform Spectrograph (SPIRE-FTS) will 
provide a spectroscopic view of a 2-arcmin-square field at 
all frequencies between 460 GHz and 1.5 THz simultane- 
ously. The 5a detection threshold of a line in a 1-hr inte- 
gration is expected to be about 1.5 x 10~^®Wm~^. The 
slope of the counts shown in Figs 1(d) and 3(b), which are 
based on the estimated, and very uncertain, counts of [Cii] 
lines, indicate that the most efficient survey depth is greater 
than 10"^* Wm-2 for HIFI and about 5 x 10"^^ Wm'^ for 



SPIRE-FTS. At a depth of 5 x 10"^*^ Wm"^ the optimal 
detection rates using HIFI and SPIRE-FTS are expected to 
be about 6 x 10~^ and 7 x 10^'' hr^^ respectively, and so 
the detection rate of line-emitting galaxies using FIRST is 
expected to be quite low. 

The currently quoted sensitivities of FIRST-illFI are 
about seven times better than those of the heterodyne in- 
struments attached to the 2.5-m SOFIA airborne telescope 
(BeckHn 1997; Davidson 2000). Thus, when the larger field of 
view of SOFIA is taken into account, a line-emitting galaxy 
could be detected using SOFIA about every 4500 hr. Hence, 
SOFIA could not carry out a successful blank-field line emis- 
sion survey. Note, however, that the instruments aboard 
SOFIA will be upgraded throughout its 20-ycar lifetime, and 
so the development of multi-element detectors and innova- 
tive wide-band spectrographs may change this situation. 

The proposed space-borne far-infrared interferometer 
SPECS (Mather et al. 1998) has very great potential 
for resolving and obtaining very high signal-to-noise spec- 
tra of high-redshift galaxies. Although predictions of the 
detectability of fine-structure lines at frequencies greater 
than 600 GHz are very uncertain, with a field of view of 
about 0.25 deg^^ at 650 GHz and a 5a sensitivity of about 
m in a 24-hr integration, about 5 [Cll] galaxies 
could be detected per day using SPECS in a 8-GHz-wide 
band (Fig. Id). Hence, although SPECS is predominantly 
an instrument to study known galaxies in great detail, its 
sensitivity is sufficient to carry out successful blank-field line 
surveys. 



4.5 Summary of prospects for line surveys 

The most promising instruments for future surveys of CO 
and fine-structure lines from high-redshift galaxies are the 
ALMA interferometer array in the 140- and 230-GHz bands, 
which should be able to detect about 15 galaxies an hour in 
a survey to a line fiux density of lO"'^*' Wm"'^ in a 16-GHz- 
wide band. Other instruments for which detection rates are 
expected to exceed one source per hundred hours are the 
future SPECS space-borne interferometer, which is likely 
to detect of order of five fine structure line-emitting galax- 
ies per (lay, and a ground-based, wide-band, background- 
limited millimetre-wave grating spectrograph, which should 
be able to detect a line-emitting galaxy every 20 hr, the same 
rate as the upgraded K- and Q-band VLA. 



5 CONCLUSIONS 

We draw the following conclusions 

(i) We have predicted the surface densities of high- 
redshift galaxies that emit molecular rotation and atomic 
fine-structure line radiation that is redshifted into the mil- 
limetrc/submillimctrc waveband. The results depend on 
both the excitation conditions in the ISM of the galaxies 
and the evolution of the properties of distant gas-rich galax- 
ies. We incorporate the latest results of CO-linc observa- 
tions and the redshift distribution of galaxies discovered in 
submillimetre-wave continuum surveys to provide a sound 
observational basis for these predictions. 
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(ii) The predicted counts of CO line-emitting galaxies are 
probably accurate to within a factor of about 5 for lines 
with J <, 7. The uncertainties are inevitable, and caused 
by the lack of knowledge of both the cosmological model 
and the excitation state of the emitting gas in the ISM of 
high-redshift galaxies. 

(iii) The predictions of the counts of atomic fine-structure 
lines in high-redshift ultraluminous galaxies are more weakly 
constrained by observations, and we present likely order-of- 
magnitude estimates of the counts in six such lines. 

(iv) The most efficient frequency for a survey that aims to 
detect CO emission is probably in the range 200-400 GHz, 
which includes the 230- and 350- GHz atmospheric windows. 
The optimal frequencies for a fine-structure line survey are 
probably in the range 400-800 GHz, which includes the 670- 
and 850-GHz atmospheric windows. 

(v) There are excellent prospects for using a range of 
millimetre- and submillimetre-wave instruments that are 
currently under development to conduct blank-field surveys 
for the redshifted emission lines from the ISM in high- 
redshift galaxies. ALMA will probably prove the most capa- 
ble facility for blank-field surveys, and should detect up to 
about 15 line-emitting galaxies per hour. The most efficient 
survey is likely to be made in the 140- and 230-GHz bands 
to a an integrated line fiux density of a few 10"^" Wm~^. 
The required sensitivities could also be achieved by future 
wide-band spectrometers on single-antenna telescopes. Low- 
J CO transitions will be detectable using the upgraded VLA 
at longer millimetre and centimetre wavelengths. 

(vi) A wide range of millimetre/submillimetre-wave lines 
should be detected in the spectra of galaxies selected in other 
wavebands, and in blank-field dust continuum surveys. Red- 
shifts for millimetre/submillimetre-selected galaxies could 
be determined directly in these wavebands if a wide band- 
width of order 100/(1 -I- z) GHz was available, with a centre 
frequency of about 200 GHz. 
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